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Introduction 


The purpose of this study Is to examine the potential benefits and 
costs of optimizing both the structural stiffness and the active control 
of aircraft in a rational manner. The ultimate goal of this effort Is 
to arrive at a unified treatment of structural and active control design 
for the stability augmentation of flexible aircraft. 

Three separate efforts have taken place during the past six months 
of effort. The first effort is an exhaustive literature evaluation in 
the area of passive tailoring for aircraft performance. During this 
effort, several valuable and previously unrecognized tailoring studies 
were uncovered. This survey was combined with similar work by Mssrs. 

M.H. Shirk and T.J. Hertz of the Air Force Wright Aeronautical Labora- 
tories to produce a paper presented at the 25th AIAA Structures, Struc- 
tural Dynamics and Materials Conference in Palm Springs, California in 
May 1984. 

The second effort involved the identification of a mathematical 
technique to be used for aeroservoelastic tailoring studies. A promising 
candidate method has been identified and is described in the following 
section. . 

Finally, two analytical models, one elementary, the other sophisti- 
cated, have been developed to illustrate the potential for aeroservo- 
elastic tailoring. Both models have essential features of "real-world" 
hardware, yet the physical understanding is not buried in a myriad of 
detail. These models are also described in the next section. 
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The Use of Structural Gains as Design Parameters 


There are two major obstacles to simultaneous treatment of the struc 
tural stiffness design optimization problem and the active controls 
problem. The first difficulty arises because of the dissimilarity of 
design variables in the two problems. This difficulty has been overcome, 
at least at the elementary level, by the selection of a characteristic 
set of nondimensional parameters for beam-like and plate-like structures. 
The state space model of an aeroelastic system can be written as: 


x = Ax + Bu 

(1) 

y = Cx 

(2) 


with x as the n-dimensional state vector, u is an m-dimensional control 
vector and y is the output vector, while A,B and C are constant coeffi- 
cient matrices. If a linear, full-state feedback control law exists, of 
the form, 

u = -Gx (3) 

then the modified system equations are: 

x = (A-BG)x (4) 

On the other hand, the equations for a structural system with passive 
control may be written as: 

x = Ax - ipAx (5) 

where is a nondimensional parameter related to stiffness cross-coupling 
provided by structural tailoring and A is a modification to the A matrix 
provided by changes in the stiffness matrix. 
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Equation 5 resembles Eqn. 6 In that 


*Ax ' BGx (6) 

If there are several variables, corresponding to tailored bays of a 
wing for Instance, Eqn. 6 becomes 

j^.A.x = BG X (7) 

Theoretically, one should be able to construct a structural modification 
in terms of ^A^ to furnish the same equivalent (in terms of eigenvalues) 
system as the actively controlled system. A major problem arises, however, 
because the elements of A are not free parameters while the elements of 
G are. Thus, standard optimal control procedures (for instance, pole 
placement) do not have an obvious adaptation. Attempts at such adaptations 
over the past six months have not proved productive. 

Fortunately, a method developed by Newson and Gilbert offers at 
least a preliminary approach to the simultaneous design problem. If 
the cross-coupling parameter \p is treated as a design parameter that is 
held fixed during the control design, its effect on the control system 
performance can be assessed by employing optimal sensitivity techniques. 
With this technique, a cost functional, J, is minimized to obtain the 
"optimal" control law for the system. The parameter is then treated 
as a design variable so that the change in 0 with respect to ^ can be 
computed usirg an adaptation of the Newson/Gilbert approach. This 
adaptation is described in the Appendix to this report. In addition, 
the sensitivity of other aspects of the control law design to 41 may 
be assessed. 
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The "important", or at least new, aspect of this development is that 
a technique is available to assess the impact of structural tailoring 
upon active control design. The models on which this study will focus 
are described in the next section. 
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An Elementary Model for Aeroservoelastic Optimization 




To study the problem of aeroservoelastic optimization, one may begin 
either at an advanced or an elementary level. The operative term would be 
"state-of-the-art." Reviews of the literature and past experience have 
convinced us that a first look at aeroservoelastic optimization (ASEO) should 
begin with an example that is simplistic, but meaningful. The model chosen 
is shown in Figure 1. This model consists of a typical section free to 
pitch and plunge as a rigid body. The design variable is, in this case, 
equal to e/b. The dimension e/b measures the distance between the static 
aerodynamic center at the quarter-chord and the plunge spring position on 
the airfoil. For a fixed ratio R = K Q /K h and with the airfoil c.g. position 
fixed, the divergence speed of the fixed root airfoil declines with increas- 
ing e/b. On the other hand, the flutter speed increases with increasing 
e/b. This provides a design trade-off situation for which an optimum value 
of e/b exists to maximize the aeroelastic stability of the system. 

If the airfoil is attached to a fuselage element that is, in turn, free 
to pitch and plunge, the situation becomes more interesting because the 
value of e now determines the attitude stability of the aircraft and values 
of e that maximize the stability of the fuselage/wing combination may differ 
significantly from those which were found for the wing alone. 

The addition of the control surface to the model provides additional 
design options. With R fixed the control effectiveness is unchanged by 
changes in e/b. Thus any design benefit or degradation is unrelated to 
control effectiveness in this idealization. 

For fixed values of the system structural and inertial parameters, 
an optimal control law may be generated. A sensitivity analysis will then 
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Figure 1 2-D ASEO Model 







be performed to assess the effect of a change In e/b (6e/b) on the active 
control of the vehicle. From this information, a new value of e will U 
selected, together with new control parameters. One possible limiting 
case of this procedure is that the active control could disappear entirely, 
meaning that passive control is sufficient to handle the stability problem. 

This model has most of the structural dynamic characteristics of an 
actual vehicle. The potential for strong rigid body/wing interaction 
exists, as does \.he capability of studying the differences between control- 
ling the stability of the wing itself (in a fixed fuselage condition) and 
the wing/fuselage combination. The most serious limitation of this model 
is the limited number of degrees of freedom. 

This analytical model does have advantages. It is a valuable learning 
tool, uncluttered by a myriad of numbers. Each step of the A5E0 procedure 
is easily understood and interpreted in light of the substantial amount 
of Information available on 2-D sections. 

Because of the limitations of the 2-D model, the AiEO procedure will 
be extended to a realistic wing with a structurally tailored span. In 
this case the cross-coupling parameter ^ is the design variable. As 
before, the objective will be to improve overall performance in light of 
lessons learned with the 2-D model. 
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Appendix 


Optimal Control and Sensitivity Derivatives 
for the Redesign Problem 

The aeroelastlc equations of motion may be written 
X c AX + Bu . 

If the control Is a linear, measurement feedback control, then 

u = GMX 

where, M = state measurement matrix (i.e. z = MX) 
and G = feedback gain matrix. 

Then the control -augmented system matrix is 

A + = A + BGM, so 
X = A + X . 

The subscript "plus" sign denotes augmentation. 

A quadratic cost function, used in linear regulator design 

( 21 ) 



J «= f [X*C*QCX + u*Ru]dt 
■'o 

where. 

C = output matrix (I.e. y = CX), 


Q * output weighting matrix, 

and 

R * control weighting matrix . 


(B.l) 

(B.2) 


(B.3a) 
(6 3b) 

is 

(B.4) 


The well known solution for the optimal control that minimizes J, 
subject to the constraint, eqn. B.l, is 


u - -R’VpX 


(B.5&) 


or, GM «= -R"*B*P (B.5b) 

where P Is the solution to the steady-state matrix Rlccatl equation, 

PA + A*P + C*QC - PBR _1 B*P - 0 (B.6) 


For computation of sensitivity derivatives, It Is assumed that the 
optimal control, eqn. B.5a, has been determined for a baseline configur- 
ation and that the weighting matrices used In this determination, Q and 
R, are "frozen" (I.e. insensitive to the design parameters, p^, so that 

4jj- = 0, and -§jj- =0). Also, the control input matrix, B, Is considered 
d P-j «P.j 

to be dependent upon the type and geometry of the control being used, and 
not upon the design parameters. So, = 0, also. Since the cost 

®P.j 

function, defined in eqn. B.4, is what determines the optimality of the 
control design. It will also be the measure by which subsequent 
redesigns are judged. 

First, the cost function is decomposed into its regulation and 
control parts. 


J = J + J 
x u 


Now, 


x ;Vo • 


where S x satisfies 


s x \ + AJS X ♦ C*QC « 0 , 


and 


J * X*S X 

u A oVo * 


(B.7) 


(B.8a) 


(B.8b) 

(B.9a) 




where S y satisfies 

v+ 4 A : s u 4 pbr ‘ 1b * p a 0 


(B.9b) 


X Is the Initial condition (time, t. Is zero) on the state vector, 
o 

The regulation cost sensitivity with respect to Is found by 
differentiating eqns. B.8a and B.8b so that 


ap. 


aS. 


X* 


O 3Pj 


where 


If* 

}p. 


satisfies 


(B.lOa) 


aS aS 3A 3A* ar 

d 1 ♦ AJ tjt + < S , -XT * vT Q C + C * Q !§-) ■ 0 • < B - l0b > 


Similarly, the control cost sensitivity can be found from 


3J 


as. 


ap 


x* 


1 


o ap. o 


(B. 11a) 


as. 


where ^ satisfies 
ap i 


as 


as. 


aA, 3A? 


A ♦ AI -p- + fs ^ S + ~ BR _1 B*P + PBR _1 B* ) * 0 
ap 1 ♦ + ap i • u ap 1 ap 1 u ap i -p^ 

(B. lib) 


Then, In general, any desired change In the costs can be effected within 
the theoretical limits of the parameters, p^ (and provided there are a 
sufficient number, NP, of parameters), as 


and 


(B.12b) 


aJ 


NP 3$ 

y x* t~ x &Pj 

,‘l 0 »P 1 o M 


Since only first order derivatives are being used, It would be wise If 
the Apj’s are kept small throughout the redesign Iterations. 

To complete this derivation, It Is necessary to obtain expressions 
for -rr- and found in eqns. B.lOb and B.llb, which are as yet 
undetermined. By first differentiating the Riccatl equation, eqn. B.6, 
and using the definition for A + , eqn. B.3a, and the solution for GM, 

a p 

eqn. B.5b, an equation that car* be solved for can be obtained, namely, 

ap^ 

lE_ a + A* 4 (-£— OC 4 C*0 — - 4 P — ^ ■ 4 p) . o (B ^3) 

8 P . ▼ a 4 S p. 3p 1 1 y 3P 1 ap i ap 1 ; { * ' 


Equation B.13 Is similar to one derived in (22) except that the Q, R, 

and B matrices are assumed insensitive to p^, and the output matrix, C, 

Is included explicitly. Note that -f— is known (see sections 2 and 3). 

3 P< »A 

dC ' 

It is assumed that t— Is known.* Then, t— can be found from 

3p i 3p i 


3p, 


• BR " 1b * I T" 


(B. 14) 


w jr 

The output matrix, C, can either be Insensitive (i.e. - — £ 0), or be 
. ap^ 

some other known function of the parameter. For Instance, if the output 
to be regulated, y * CX, consists of Internal structural loads, then 
C will resemble some portion of the structural stiffness matrix. 
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